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Abstract. Galaxy clusters are sources of gravitational ra- 
diation. The main aim of this paper is to give numerical es- 
timates and theoretical description of the relevant features 
of the gravitational radiation coming from an appropriate 
spatial distributions of galaxy clusters. Since no analytical 
approaches are currently available to describe the strongly 
nonlinear regime, our numerical approach -combining nu- 
merical simulations with statistical arguments - seems to 
be an useful way of studying the main features of that ra- 
diation. Although far to be detectable with present tech- 
nology, we advance some ideas about future observational 
strategies and its cosmological implications. 
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1. Introduction 

The cosmological stochastic gravitational- wave background 
produced by mildly nonlinear evolution of density fluc- 
tuations can be estimated using analytical perturbative 
methods (see, e.g., Matarrese & Mollerach 1997 and refer- 
ences cited therein). However, the estimation of that emis- 
sion, in the strongly nonlinear regime, would require non- 
perturbative approaches which are not currently available. 
Galaxy clusters are evolving beyond the mildly nonlinear 
regime and they undergo the so-called virialization process 
(violent relaxation). Then, the following question arises 
(Matarrese & Mollerach 1997): have the clusters produced 
a significant amount of gravitational radiation? Since the 
problems with the analytical calculation of the gravita- 
tional radiation coming from strongly nonlinear galaxy 
clusters has not been still solved, we have estimated that 
radiation numerically. 
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In a previous paper (Quilis et al. 1998b), the gravita- 
tional radiation released during the formation of an iso- 
lated galaxy cluster was analyzed. In this paper, we have 
extended that calculation to the case of the emission com- 
ing from spatial distributions of galaxy clusters, thus we 
have performed, for the first time, a quantitative nonlin- 
ear calculation of this emission. Hence, we have obtained 
a numerical estimate of the contribution to the stochastic 
gravitational-wave background due to a given spatial dis- 
tribution of galaxy clusters. Although the resulting emis- 
sion has appeared to be weak in our calculations, we want 
to stress that, at least from the theoretical point of view, 
it is worthy to know the main features of that radiation. 

Cluster simulations including both baryonic and dark 
matter have undergone a high level of development (e.g. 
Kang et al. 1994, Metzler & Evrard 1994, Navarro et al. 
1995, Anninos & Norman 1996, Gheller et al. 1998). The 
baryonic component is evolved either with an Eulerian or 
with a Lagrangian code, the treatment of the dark matter 
component is based on appropriate N-body techniques, 
and both components are gravitationally coupled through 
Poisson's equation, which is usually solved using the Fast 
Fourier Transform. Simulations based on a 3D Eulerian 
code (Quilis et al. 1998a) were used to estimate the gravi- 
tational radiation generated by individual galaxy clusters 
(Quilis et al. 1998b). Computations were performed in 
the framework of the standard Cold Dark Matter (CDM) 
scenario. The Hubble constant and the density param- 
eter were assumed to be Hq = 50 Kms^^Mpc^^ and 
rig = Ij respectively. The same assumptions and method 
have been used in the cluster simulations used in this pa- 
per (this procedure is justified below). Appropriate initial 
conditions have been chosen to get both rich and stan- 
dard galaxy clusters. Our rich clusters have a X-ray lumi- 
nosity (Lx) of the order of ^ 10'*'' erg/s, a temperature 
T > 3 X 10^ K and, a total mass inside the Abell radius 
(3 Mpc) M > 4 X 10^"* Mq. Our standard cluster has 
Lx - lO'^^ gj.g/s^ T < 6 X 10^ and, M < lO*-* Mq. 

As it was discussed in Quilis et al. (1998b), the gravita- 
tional radiation from a galaxy cluster produces progressive 
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deformations on some material systems, this is what was 
referred to as the secular effect. In the particular case of a 
system formed by two test particles, the deformation re- 
duces to a relative variation of their separation distance. 
Since this variation is proportional to time, then, leaving 
aside the question on the characteristic frequency (cru- 
cial for detectability), it reaches values which are in the 
range of space-based laser interferometric observatories of 
gravitational waves. 

There are various conclusions of Quilis et al. (1998b) 
to be taking into account here: (i) The secular effect is 
approximately proportional to D~^, where D is the dis- 
tance from the cluster to the observer, (ii) clusters located 
at a distance D > 600 Mpc only produce a small effect 
which is neglected in the simulations of this paper, (iii) 
for distances D < 600 Mpc, gravitational waves can be 
considered as propagating in the Minkowskian space tan- 
gent to the Friedman-Robertson- Walker spacetime at the 
emission point. Although the secular effect is only signif- 
icant in the near zone {D < 600 Mpc), no differentiation 
is made - along the paper - between this effect and grav- 
itational waves (in the radiation zone), this is because, in 
both cases, the deformations produced on detectors have 
the same sources (galaxy clusters) and the same transverse 
signature and, also, because these deformations are esti- 
mated with the same formulae holding in both the near 
and the radiation zones, (iv) the secular effect produced by 
each cluster does not appear isolated, but superimposed 
to the effect due to other clusters, and (v) the internal 
dynamics of clusters having similar features (L^, T, and 
M) is expected to depend on the cosmological parameters 
(fig, Hq .,...) weakly, while the spatial distributions of clus- 
ters is sensitive to these parameters. Consequently with 
this last point, we have fixed Qq and Hq for simulating 
individual clusters, while the chosen spatial distributions 
mimic some features of cluster catalogues. 

2. Simulation strategy and results 

Only very crude simulations of the total secular effect pro- 
duced by cluster distributions were given in Quilis et al. 
(1998b), where some numbers were found with the essen- 
tial aim of proving that the effects of many clusters do 
not cancel among them. The main goal of this paper is, 
exactly, this one: to present improved calculations show- 
ing the features of the signal coming from an appropriate 
spatial distribution of galaxy clusters. This spatial distri- 
butions of clusters could be chosen either mimicking the 
observed one or simulating them in a certain scenario of 
structure formation in the Universe. In any case, the distri- 
bution can be restricted inside a sphere of 600 Mpc radius. 
According to Dalton et al. (1994), the analysis of the APM 
galaxy survey leads to a mean density of rich clusters of 
4.25 X 10~^ Mpc~ and a two-point correlation function 
of the form (cc — {ro/rY, with tq = 28.6 Mpc. All the dis- 
tributions of rich galaxy clusters used in this paper - to 



estimate the total secular effect - have been constrained 
(Pons-Borderia et al. 1999) to have these features. 

In order to give a complete and clear physical descrip- 
tion of the secular effect, a simple detector formed by two 
test particles, A and B, is appropriate. Let us analyze the 
response of this detector to the gravitational waves from a 
cluster (secular effect). In the Transverse Traceless (TT) 
gauge, the relative motion of these particles is fully given 

TT 

(Misner et al. 1973) by the quantities h^^ , which describe 
a propagating small perturbations of the space-time struc- 
ture (gravitational wave) . Direction x'^ = z coincides with 
the line of sight of the cluster producing the gravitational 
radiation, while directions a;^ = a; and x^ = y are in a 
plane orthogonal to this line. There are four nonvanishing 
components of /i^ ■ satisfying the relations h^. 



y - hy^ and 

TT 

and h^y = h^ are 



Ky = ~'*^^' ^^"^^^^ only hZ = 

independent quantities defining two polarization states. In 
TT gauge, there is a system of coordinates attached to A, 
in which the coordinate of the particle B undergoes the 
following variations from present time ip to time io + At: 

Ax^-{xoh+{to) + yoh^{to)}A{t) (1) 

Ay ^ -{xoh^ito) ~ yoh+(to)}A{t) (2) 

Az = (3) 

where At is much smaller than the period of the gravi- 
tational waves emitted by clusters {At << 10^ yr), the 
subscript " 0" stands for the initial coordinates of the par- 
ticle B at time to, the overdot stands for a time derivative 
and, the quantities /i+ and h^ are computed at point A 
and at present time. 

Each cluster produces a secular effect described by 
Eqs. dll) - (pi). These equations show that the effect is 
proportional to At. The estimate of Ax and Ay requires 
the knowledge of /i+ and h^ for the chosen cluster. An 
explicit computation of these quantities is only possible if 
cluster evolution is known, as it occurs in the case of nu- 
merically simulated clusters. The total secular effect pro- 
duced by a distribution of clusters can be obtained from 
the /i_|_ and hx quantities corresponding to each cluster. 
Obviously, the simulation of each one of the ~ 4000 rich 
clusters located inside a sphere of 600 Mpc radius is out 
of current computational capabilities; therefore, some type 
of statistical treatment of the problem is necessary. Quan- 
tities /i+ and hx must be assigned to each one of the 
clusters belonging to some spatial distribution using ade- 
quate criteria. Let us motivate these criteria listing various 
considerations: 

(a) The gravitational radiation from a cluster is the 
superposition of the radiation produced by the motion of 
many particles of dark and baryonic matter inside the clus- 
ter; hence, low levels of polarization are expected and the 
assumption that quantities /i+ and hx are independent is 
good enough. Furthermore, clusters radiate incoherently 
and, consequently, the waves from different clusters have 



Vicent Quilis et al.: Gravitational waves from galaxy cluster distributions 



distinct uncorrelated phases. At emission time, each clus- 
ter should be in an evolution state different from and in- 
dependent on the state of any other cluster. 

(b) Several simulations of rich clusters have been an- 
alyzed in order to define an interval where the values of 
ft.+ and hx are distributed. Let us assume a sphere of 
100 Mpc radius centered in one of our simulated clus- 
ters. We could observe the central cluster from any point 
on the sphere. The radius passing through the observa- 
tion point would be our z axis and, taking into account 
the evolution law of the cluster - which has been numeri- 
cally simulated -, quantities /i+ and h^ can be computed 
for each observation point. Then, a great number of these 
points can be considered and the maximum and minimum 
of these quantities can be easily estimated. The mean of 
the rnaxima obtained from all the simulated rich clus- 
ters, h '^ {h+{inax)) ~ {hximax)), has appeared to be 
~ 9. X 10"^"'^ yr~^; hence, we assume that quantities h+ 
and hx are distributed in the interval {—h, h). 

(c) The same study has been done for standard clus- 
ters. In this case, the resulting h value is ~ 1. x 10~^^ yr~^; 
hence, rich clusters produces an effect which is one order 
of magnitude greater than that of the standard clusters. 
On account of this fact, only the distributions of rich clus- 
ters are considered in this paper; nevertheless, more work 
should be done to estimate the contribution of standard 
and small clusters, which are abundant structures produc- 
ing weak secular effects. 

After these comments, the following method seems to 
be appropriate to give values to h+ and hx'- According 
to (a), these quantities are generated as statistically in- 
dependent numbers and, on account of (a) and (b), each 
of these quantities is assumed to be a random number (ry 
for /i_|_ and ^ for hy) uniformly distributed in the interval 
(-1,1) multiplied by the mean value h. Finally, if a given 
cluster is not located at 100 Mpc from the observer, but at 
a distance D, number h must be multiplied by the factor 
lOO/I?. These criteria plus Eqs. (|l|) - (H) allow us to com- 
pute the relative variation of the distance AB produced 
by an arbitrary cluster located at distance D (in Mpc) 
from the detector. After trivial algebra, the following key 
equation is easily found: 



I 



D ^' 



where 9 is the angle formed by the segment AB and the 
cluster line of sight. For = 0, the segment AB is aligned 
with the line of sight and no deformation is produced at 
all. In the case = 7r/2, the segment AB is orthogonal to 

a; 

the line of sight and the cluster produces a maximum — — 

independent on the orientation of the AB segment in the 
plane orthogonal to the line of sight. 

Given a spatial distribution of clusters and an orien- 
tation of the segment AB, Eq. (0) allows us to assign a 

M 
small — - to each cluster and, then, all these values must be 



added to find the total secular effect produced by the clus- 
ter distribution. If this distribution is not altered, but the 
orientation of the segment AB is changed, the total secu- 
lar effect measured by the detector changes (anisotropy) . 
According to Eq. (0), this change occurs because numbers 

77, ^, and h keep unaltered, but the angle corresponding 
to each cluster varies. Given a spatial distribution of rich 
clusters, we have all the ingredients necessary to generate 
full maps of the sky; namely, maps where all the AB ori- 
entations are considered. These maps completely describe 
the anisotropy of the total secular effect. 

One map of the full sky is displayed in Fig. 1, where 
quantity /i = '^^St ^^ given for all the directions join- 
ing the center of a sphere (A) with the points (B) lo- 
cated on a hemisphere. There are anisotropics in the sense 
that, for the chosen cluster distribution, the /i value - 
displayed by the grey scale- depends on the orientation 
of the segment AB in the space. The map of this Figure 
has been expanded in spherical harmonics turning out in 
the superposition of a monopole ao,o = —8.17 x 10~^° 
and the following independent quadrupole components: 
the real component 02.0 ~ —1.38 x 10"^", and the com- 
plex ones a2,2 = a^ -2 = (1-44 x 10-20,-1.78 x lO^^O) 
and 02,1 = -a^-i = (-4.25 x 10-20,1.81 x lO-^O). Any 
other multipole appears to be negligible. These numeri- 
cal results are easily understood taking into account Eq. 
(P, where the dependence on sin^9 shows that the effect 
of each cluster is the superposition of a monopole and a 
quadrupole. We can conclude that the total secular effect 
can be completely described by the six components of the 
quantities (ao,Oj ^2,2, 0^2, 1, 02,0)1 which depend on the cho- 
sen cluster distribution. 

Various simulations have been done using different dis- 
tributions of clusters. Many orientations of the segment 
AB have been considered, and the mean and standard de- 
viation a of the predicted /i values have been calculated. 
The mean value of fi changes from simulation to simula- 
tion and the standard deviation is a rather stable quantity. 
This quantity, which measures typical deviations with re- 
spect to the mean and, consequently, anisotropy, ranges in 
the interval (1. x 10-^° yr-\ 2. x 10-^° yr-^). In the worst 
case, 4^ = ^At takes on the values '^ 10-^^, ^ 10-^^ and 
~ 10-^^ in time intervals of 30 s, 3.6 days and a decade, 
respectively. 



3. Discussion and speculations on observational 
strategies 

The efficiency of the gravitational wave emission from a 
rich cluster is very low, « 10"^^ (see Quilis et al. 1998b). 
This means that during all the age of the Universe a single 
cluster would radiate a gravitational energy of « IQ-^^ 
times the cluster mass. The contribution of that energy 
to the present density parameter, fig^,, would be of the 
order of w 10"^^ flc, being ilc the contribution of all the 



Vicent Quilis et al.: Gravitational waves from galaxy cluster distributions 




-5.53 



2.23 



|Ll 



Fig. 1. Map of the secular effect pro- 
duced by a distribution of rich clus- 
ters. The quantity /i x 10'^", with 
H = l~^-^ in units of yr""'^, is plot- 
ted for a semisphere. Lambert's pro- 
jection is used. 



clusters to the density parameter. That rough argument 
suggests that Vtgy^ would range between 10"^"^ - 10"-^*. A 
very low value. 

The gravitational waves from galaxy clusters produce 
anisotropics on the Cosmic Microwave Background (CMB). 
Other gravitational waves generated during inflation or 
in another early process would lead to "primary" CMB 
anisotropics because they were present at recombination 
time; however, the gravitational waves coming from galaxy 
clusters were emitted after recombination and they can 
only produce "secondary" anisotropics, which are due to 
the motion of the CMB photons in the time varying grav- 
itational field associated to these waves. Although these 
anisotropics exist (detailed estimates are in progress) , they 
are expected to be too small -for detection- due to the low 
values of ^gw given above. This means that observations 
of the CMB seem not to be appropriate for detecting the 
gravitational background from galaxy clusters. Could we 
use interferometry, as in standard detectors of gravita- 
tional waves, to detect this background? 



3.1. Some clues for future observational strategies 

Measurements of the total secular effect do not require a 
fixed orientation of the AB segment; suppose, for instance, 
that particles A and B move with the terrestrial equator; 
then, the AB direction depends on time in a well known 
way and it covers all the directions of the equatorial plane 
during a day. Every day, the total relative variation A/// 
is given by the integral 



(A//Oday=-^y'%(/3)d/3 



(5) 



where angle /? defines the AB direction inside the equa- 
torial plane and W is Earth's angular velocity. We have 
fixed a cluster distribution and, then, various planes play- 
ing the role of the equatorial one have been considered. 
See, for instance. Fig. 2, where functions n{0) are dis- 
played for three of these planes. Using a large number 
of planes and Eq. d|), many possible values of (AZ//)day 
have been calculated. The mean of these values and the 
standard deviation with respect to it are 4.93 x 10"^'^ and 
1.8 X 10^^^. Other cluster distributions have been consid- 
ered and the means and standard deviations appear to be 
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Fig. 2. Plot of the quantity /ix 10'^° (in yr^^) as a function 
of the angle /3 (in degrees) defining the orientation of the 
AB segment inside a fixed plane. Three different planes 
have been chosen at random. 

different; nevertheless, the orders of magnitude keep unal- 
tered. We can conclude that the detector would undergo 
successive pulses with amplitudes of the order of 10^^^ 
and periods of 6 hours (the period is T = tt/2W, see Fig. 
2). The pulse appears as a result of the detector rotation in 
the anisotropic gravitational field created by the clusters. 
Since the sensitivity of laser interferometry-based de- 
tectors decreases with the signal frequency v, they can- 
not observe -directly- secular effects with almost negligi- 
ble v ; nevertheless, as discussed above, rotating detectors 
would receive effective pulses - with frequency v — 2W/t: 
— whose possible detectability would require further anal- 
ysis. The angular velocity of the detector fixes the pulse 
frequency 2W/'k and its amplitude (which is proportional 
to W~^ according to Eq. (||)). According to our estimates, 
the amplitude reaches the order 10^^^ for pulse frequen- 
cies of the order of 10~^. As the angular velocity increases, 
the pulse frequency increases - gaining in sensitivity - 
but the pulse amplitude decreases. Taking into account 
these considerations and the fact that interferometry is a 
very accurate technique, the design of specific experiments 
for detecting the above systematic and repetitive pulses - 
produced by the clusters on rotating detectors- deserves 
attention. 



above) . Distinct values of these parameters would lead to 
different spatial distributions of clusters and, then, to dis- 
tinct secular signals. For instance, the dependence of the 
secular effect on the Hubble constant can be easily ana- 
lyzed. As this constant changes, all the distances are mul- 
tiplied by the factor Hq^, and the angles (9) between the 
AB segment and the cluster lines of sight do not change; 
hence, the secular effect, which is roughly proportional to 
-D~^, is approximately proportional to Hq. The secular 
signals could be predicted for a wide range of values of 
the involved parameters and, then, comparisons of the re- 
sulting predictions and future observations could give new 
bounds on some of these parameters. These bounds would 
be very interesting from the cosmological point of view; 
however, taking into account that: (i) the secular effect 
is small and, (ii) its detection is not expected to be easy, 
our considerations about the cosmological consequences 
of detecting the gravitational radiation from clusters have 
to be read wisely. Other observations could be much more 
effective to constraint the cosmological parameters; for in- 
stance, CMB observations from future spatial missions as 
MAP and PLANCK. 
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3.2. Cosmological consequences 

The secular effect and its anisotropy are expected to be de- 
pendent on the value of some cosmological parameters (see 



